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ABSTRACT 


A nonlinear clutch type servomechanism is proposed to position 
ivescontrol surfaces of an underwater vehicle in accordance with specified 
merformance criteria, 


Two system configurations are discussed and compared; 


(1) A two-state system employing pulse width angular velocity 
Medulation On Counter rotating clutches. 


(2) A three-state system employing counter rotating clutches and 
a brake with error sampling. 


Stcady state load vibration will be present in each case. The 
degrading effects of these vibrations are developed and discussed. 


Clutch power dissipation phenomena are discussed. 


Actuator specifications are examined in light of overall system 
dynamics. 


The present state of clutch technology is examined to fulfill system 
Gecduirements. 


iitesis Supervisor: George ©. Newton, Jr. 
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CHAP TE Rt 


INTRODUCTION AND SUMMARY 


ieee ROBE MM STATEMENT 

Switching amplifiers, in general, offer higher efficiency and 
higher power per unit weight and volume than do continuous amplifiers. 
This research investigated the application of clutch type actuators 
Swetated in a switching mode to the feedback position control of the 
control surfaces of an underwater vehicle. The vehicle was assumed 


Eenbe Operating at maximum speed. 


1,1 SYSTEM CONFIGURATION 

A self-propelled underwater vehicle when properly designed with 
neutral stability in yaw-rate will exhibit fastest response in that mode 
of motion. A control surface position servo satisfying the requirements 
of yaw-rate tracking will be adequate for roll and depth control. 

A linear small signal model for the yaw-rate loop is indicated 
ioFic, 1, l(a). 


fice joadividual@transter functions are 


Yo _A(sta) 
6 ~ s(stb)(stc) 
§ ne 
ae gee ee ee 
a a + 20w s ron 
n n 
where Gear 54 are preassigned gain factors 


ge ire 


Wa emmaie SUCcGmependent yvenicle parameters 
O, = natural frequency of servo actuator system 


ee 
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Yi + oF Rudder 5, Vehicle % 
€ . =| Servo “| Dynamics > 
7 Rate [L 
Gyro 
Fig. 1.1 (a) Yaw-Rate Loop 
id —— So 
ike —_ 
Ks = TORQUE AMPLIFIER GAIN FACTOR 
: 2 
in -FF sec 
Jy = 0.08 ee LOAD INERTIA WITH ADDED 
HY DRODYNAMIC ENTRAINMENT 
ins sec 
By = 5 ee LOAD VISCOUS AND HYDRODYNAMIC FRICTION 
in +f | 
KE = 1910 ah LOAD RESTORING TORQUE SINCE CONTROL 
SURFACE SUPPORT IS FORWARD OF QUARTER CHORD 
Ke = RATE FEEDBACK GAIN FACTOR TO GIVE 0.5 DAMPING RATIO 


Fig. 1.1(b) Control Surface Actuator 





oe 


ux 
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damping ratio 


—- 
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command yaw rate 


oF = COmiimend ruddaer detlection 


bo vehicle yaw rate 
ae actual rudder deflection 
A linear model for the control surface position control system is 
indicated in Fig. 1.1(b). The control surface system frequency response 
has been specified as flat to 7 cps with -3dB point at 27 cps. (4) a ive 


Pomawenamplitier K. must deliver up to 1000 in -lb torque ata 


6 
maximum speed of 100 RPM. 
ie AGU ATOR CONFIGURAZGION 

The design of a linear control system for control surface positioning 
is straightforward and can be realized with an electric or hydraulic 
actuator. These devices have inherent disadvantages, notably energy 
“Oiviemolommin the former and shelf life in the latter. 

Clutch type actuators offer superior performance in many respects 
although analysis is difficult due to nonlinear describing equations, 
Counter rotating clutch actuators can be operated in a push-pull fashion 
or ina switching mode. The former is undesirable from efficiency 
and heat sinking considerations. Switching operation is considered in 
this research. 

Two basic configurations are considered: A two-state pulse width 
modulated counter-rotating clutch arrangement and a three-state 
System utilizing two counter-rotating clutches and a brake with state 


determined by sampled errors. 





aac 

The two-state system is functionally diagrammed in Fig. 1. 2(a). 

The pulse width modulator output is designed to operate ata 
Constant Geeurrence rate witheuxed output magmitude positive or negative 
maurdali cycle duration determined by input level. ) The pulse width 
modulator output drives the actuation mechanism of the clutch. 

High inertia clutch input members are driven at constant speed 
imemmriesvellcic Prime mover, Assunre that input and output members 
become instantly locked on signal from the pulse width modulator. 

The clutch output member velocity will be a square wave. The spring 
K; transforms this velocity square wave into a triangular torque wave 
Pavesce smOrt-time mean levelis transmitted to the load. The load and 
filter inertia attenuate torque components at switching frequency 

and higher harmonics. 

Compensation and feedback networks insure reset action and 
elabilize the system, Reduction pearing is necessary to reduce clutch 
51ze. 

In the steady state the clutch output velocity will be a symmetrical 
square wave. Some wobble will be inherent in the control surface load 
at switching frequency and its harmonics. 

dhe three state system operation is self explanatory in Fig. 1.2(b). 
System operation is quite similar to that of the two-state configuration. 
fine notable exceptions are the more complicated decision box and lack 
of steady state output oscillations at the expense of steady state errors, e. 
Limit cycle oscillations can occur at sampling frequency if the forward 


loop gain is high enough. 





=o 

Detailed design of the clutch actuator system is the subject of 
a separate thesis and will not be dealt with here. 

The advantages of clutch type actuators derive from the high 
torque squared-to-inertia ratios and low control power requirements. 
Smelt lite ys no problem with clutch servos. 

The disadvantages of both configurations include heat dissipation 
at contact surfaces and output wobble with the possibility of exciting 
structural resonance, 

Selection of the optimum switching of sampling frequency will 
involve tradeoffs between output wobble amplitude, which will decrease 
with switching frequency, and power dissipation in the clutches, which 
will increase with switching frequency. The fundamental switching 
frequency must lie above load resonance and away from see 


BeooOnlances, 


mo CONTROL SURFACE 

The adverse effects of control surface wobble were investigated 
and acceptability criteria were assigned, 

Control surface wobble manifests itself as a perturbation on 
the control surface angle of attack. The hydrodynamics of the situation 
feeo CXamined, Since mean incidence rémains unchanged the dynamic 
Meteot the Control surface is not affected. Viscous drag effects were 
mecoma order in magnitude, tn fact, with the proper phasing and correct 
relationship between frequency, vehicle velocity and control surface 
chord length, propulsion will result. Cavitation inception was considered 
unlikely since the disturbance velocities were insignificant compared 


with quasi-steady levels. 
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Sound generation by the oscillating control surface proved 


to be a significant result, The sonic intensity level can be 


determined from the equation 


el A? 2 
(—— ) cos @ 
dt Z 
Ber)). = — aaa 2 watts /m 
is ew fey te! r 
Oo O 


where (fe iit toOnecese~erted by mydr1oro) | 


= angle between direction of observation 
and littrorce 


pees density of undisturbed fluid 
ea speed of sound in undisturbed fluid 
r = distance from source to observer 


This relationship is developed in Chapter Two. Thus, knowing 
lift as a function of angle of attack and angle of attack as a function 
Bimtiare, theysonic intensity can be obtained. The sound pressure level (SPL} 
eombe Obtaimed trom intensity. A convenient measure of sound is the 
Bource level (far field radiation extrapolated to r= 1M) expressed in 
giecibe)] ratio to 0,002 microbar (2 x 107° Newton tivencnatn 

For a given maximum control surface rate specification (100 RPM 
in this case) the source SPL can be determined as a function of switching 
frequency in the two-state configuration, A limit cycle amplitude 
determination must be made for the three-state configuration. 


Mie two-state source Sil, was ebserved toidecrease 20 dB per 


mecade of switching frequency. 





i 
I+ Gime eco POWER DISSIPATION 

Retetrring te Fig. 1. Zia) if K, were zero,the amount of energy 
meGessaiy toraccelenate the clutch output member from full clockwise 


angular velocity is 


F) 
" 


a 
sl J & joules 


ey 
1! 


clutch output inertia 


~ 
t" 


angular velocity of input members 


Since this operation is performed twice ina switching period, 


Pe -=PZ2F of f watts 
c sw 
io =) sWitecningsireg (cps) 


inherdiscipatvon rating,ef the clutch must he abeve this level, 
The situation can be alleviated somewhat by making Ky as stiff as 
possible without exceeding clutch pullout torque in spring windup. An 
analog simulation bears this out (see Fig. 4. l(b) ). This result is 
consistent with two-state system design objectives in that clutch 
Purlout torque miust be higher than the greatest demanded. 


Epetnevamalor simulationmtne filter inertia J. (Fig. 1.2(a)) was 


f 


Bomsiaered intimte, in the actual system J, > oS. 1Omeectimcie nt 


fi 
filtering so this assumption was not a bad one. The simulation is 
metailed in Ghapter Four. 

Since power dissipation is a linear function of switching 
frequency and a quadratic function of input member speed the gear ratio 
siould be selected to minimize input member speed. Unfortunately 


Merctmuorduesrecqwirements determane the pear ratio since the clutch 


dimensions must be kept small. Pullout torque is determined by clutch size, 
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[Ptheseontro! surfac® can be pivoted nearer to its center of 
pressure the output torque required (1000 in-lbs) will be reduced 
amastne meecessary reduction ratio Can be reduced permitting lower 
clutch input member speed while maintaining maximum control surface 


Bositien rate at the specified 100 RPM. 


mo OvmRAI I SYSTEM DYNAMICS 

Since the control surface servo is embedded within the overall 
vehicle dynamic control loop, a study wascarried out to determine 
mie eCitect of actuator Seerciuilelts reduction on the dynamics of the 
total system. 

The vehicle exhibits fastest behavior in the yaw rate loop (Fig. 1. 1{(a)). 
The control surface servo system was modelled as a second order 


system with damping ratio equal to es : 


2 
2 
6 w 
BeeC ae = 
oF Sirens stew 
n n 


It can be shown that damping ratio = 0.5 minimizes integral square 
earor for a linear oer iore servo for a — inpUten) Enis behav donee am 
be achieved by proper selection of forward loop gain and output rate 
feedback. 

Bode plots of the open yaw-rate loop were constructed for the 
specified servo and three slower servos. These data were transferred 


toa Nichols chart for comparison (Fig. 1.5). 
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Fig. 1.5 Nichols Chart of Yaw-Rate Loop Frequency Response 
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It was concluded that the response speed of the control surface 
actuator could be halved without degrading the closed yaw-rate 


loop performance. 


imeeee Ud Ch SPECIFICATIONS 

Mnesstate or clutch technolopy was Scanned and unfortunately 
Pacsent procuction line clutches do not fulfill system requirements. 
Syuavenme rates in the range 75 to 150 cps are meccssany for the two 
PibetemoystcCia = Jhis gives rise to am input-output rise time of one or 
two milliseconds. This demand has been met by very small instrument 
magnetic particle clutches. Figure 1.6 compares characteristics of three 


production devices with proper torque and speed characteristics. 


Maximum Maximum Maximum Rise ee Dissipation 
Clutch Torque Speed Dimension Time > 4 Rating 
Type (in-lb) (RPM) (in) (im sec) (in-lb sec x10 ) (watts) 
solenoid- 
ectated 6g 7500 3.7 30 6.93 136 
Mech. 
Briction 
Magnetic | 
Particle 5 3400 orl 30 Lot 75 
Selenoid - 
Actuated 
Vacuum- 60 1800 5 8 9.3 Wag 
Operated 
Mechanical 


Fig. 1.6 Comparative Clutch Characteristics 





a wae 

In addition, the magnetic particle clutch exhibits 0.6 in-lb drag 
Poe dUcealiezerGec<citatlon andtme vVacuum-operated clutch requires a 
small vacuum pump. 

It may be concluded that the two-state actuator configuration is 
Hovechmredilyveteasiple at present, Since the three-state system has 
not been fully investigated, a working model may be constructed using 
Hjespresent clutch technology. High speed magnetic tape drives have 
indicated a need for fast mechanical clutches and current research may 
Bolve a device capable of meeting the requirements of the two-state 


mewuator System, 





GHAPTER il 


SOUND GENERATION 


The intensity of sound resulting from a lifting body in fluid flow 
has been investigated, we Rayleigh has shown that any frequencies in 
aeibuid tloweare equal to those of the sound generated. 

Lighthill has considered the equations of motion of fluids and 
compared them with the equations of sound propagation in a medium 


erect nhe G@xact Cquatmons ton the former may be written in the form 


oS o> wate Ts, (2.1) 
where tee = es . a ae 
1j iy 9 i) Oueein 
p = density 
Pe SCOMipGessiVessuress, Leng or 
C. es Omle velocriyetm tld at rest 


< 
t! 


VeLoc PyEcomMiponmemtin ir ection x. (ole Ces) 


Physically, this means that sound is generated by a fluid flow 
exactly as in a uniform medium at rest ani LS) Cle Gl Ony Divaercte Tina lily; 
applied fluctuating stresses. If we regard the flow as known it is possible 
to evaluate the sound generated, 

Equation (2.1) is the non-homogeneous wave equation. The most 


general solution (Stratton 1941) is 











2 
= @ ia. 
p(x) - p Sas oe dy 
4nC = ies) || - y | 
] oe: ] or ] or op 
ae r on + 2 Qn? * ©F On Ot re 
S 





= Ae 





Z 
Oo | 3 9 9 
In this equation —- J ae =O ; oP are taken at retarded times 
OY; 0Y; n ot 
= ~ fo oeGe f= |x - y | and n is the outward normal from the 


O 


fluid. The first integralis taken over the total volume V_ external 
to the solid boundaries and the second integral is taken over the 
BUriacemo Of the solid boundaries. 

Curle has shown that the solution may be simplified under the 
assumption that: 

(al) ee) a ee Geile ld. ae Ox lina tion } 

(2) the solid boundaries are rigid, i.e., non-dilatational 


(3) the flow is incompressible (Mach No. << 1) 


Under these assumptions the solution reduces to 








x 
A ] i re) 
ep FO Sy —— F(t) 
O dae xe ot 1 
O 
where F(t) = J iP) (ee t) d Sty); the total resultant force 
Ss: 


exerted on the fluid by the solid boundaries. 

In acoustic media, the sound pressure field is proportional to 
the density fluctuation field by the constant of thus 
] iO 


A 
Ie 4 mc Z ot FAC) 
ae 








This is the radiation field of a dipole source with axis parallel 
Pomr(t) the net resultant force exerted on the fluid by the boundaries. 


Rewriting in more convenient variables: 








A ] cos §@ 
p 


- 4 ac r 
O 





F(t) (2.2) 


© 
=a 


where @ is the angular distance between the radius vector r and the 
axis of F’, (t). 
Since intensity is defined as the acoustic power transmitted through 


Amumit area 


Z 
I= —= 
(2 
Povo 
wy . Z 
p = Pincatesdtatre macouscic pressure (Newtons/M ) 
PP Ors C= characteristic impedance of acoustic medium 


(MKS RAYLS) 


Generally it is more physically meaningful to consider a decibel 
natio referred to a datum level for intensity and pressure measurements. 
A customary quantitative measurement in underwater acoustics is 


the source level sound pressure level. This is defined as 


ae 
SPL = 101log,, —h;- dB 
O 


where a is far-field acoustic mean square pressure extrapolated 

to a radia] distance of one meter from the source, the reference pressure 

Beime 0.0002 microbar (2 x 10° °N/M“} fOorsuMmcdeTWabehenolseinieas Ubenne nt an 
if the force exerted by thejbody on the fluid can be determined then 


the SPL can be calculated, 





SA Pie iil 


Por Oi) yY NAMICS aime CONTROL SURE ACE OSCILLATION 


Principles of incompressible, non-cavitating hydrodynamics 
Camebe utilized to determume the external forces acting upon the 
fluid which will yield the sound pressure field. Sound generated by 
the oscillating control surface only is considered. This sound pressure 
will be additive to any presently inherent in vehicle operation exclusive 
Siether control surface system. 

imeduation (2. 2), F(t) must be determined to find the SPL. This 
Fetaistorce exerted by the body on the fluid is due primarily to dynamic 
itt#eecdirected norinal to the approach flow. A secondary force is due to 
viscous drag which is oriented peel to the approach flow and generally 
equal to about one twentieth of the lift in magnitude. 

For an oscillating hydrofoil the oscillating component of lift is a 
function of approach velocity, hydrofoil sear fluid density, and the 
magnitude and frequency of the oscillation. For a steady hydrofoil the 
Gymainic lift is a function only of approach velocity, foil geometry, 
fluid density and angle of incidence with approach flow. A frequently 
used dimensionless parameter is the lift coefficient Cre 


G S L = Ka (gn a< in) 


L e Z 
Zoo 


for a symmetrical thickness form. 


ma ae 





ae ey 


where i = . dynaimie ant 
0 = fluid density 
© = area of foil in the span-chord plane 
US se approach tilew yeloeiry 
a = radian angle of attack 
K = proportionality constant dependent on aspect ratio 


ihe aspect ratio is the ratio of span squared to foil area S. For two 
smensiamal flow {inftimte aspect ratio), K = 27. For zero aspect ratio 
(chord length much greater than span), K = 1/2. It is known that at 
mecuced frequencies preater than 0,5 the oscillating component of 


the flow becomes two-dimensional. 


K = rab Swereaucea trequcncy 
w = radian frequency of oscillation 
5 = semi-chord length 
Vo = eaperoach flow velocity 


In the course of a mathematical investigation of wing flutter 


Baenomena, nee don oem) has shown that 


(2) 
oe ea) ae Se 
stay Bie (k) + fee (k) 


fen a lifting surface in rotary oscillation. 


where cna = Complex value of oscillatory lft component 
Log = quasi-steady lift corresponding to 
Y maximum inclination 


H') = Hankel function of second kind of order n 





+~).oe 

The function C(k) is plottedin Fig. 3.1. Its behavior is quite 
cHoniceueromtiat Ot a lat metwores Phe pwedueced frequencies of the 
VenlGle"comrnro)! Surface in the range of interest were found to be greater 
than 1.0so that C(k) was assumed equal to [ 0. 35) sp yee 0] throughout 
the analysis. 

iiemonesolne Conclusions Nave been predicated on the assumption 
that the foil is oscillating in a sinusoidal] fashion with time. 

Petemiunc LO Mic lad assume thatethemclutch issperlectin that 


imesoutout member locks on with the appropriate input member at the 


switching instant. The angular velocity of the clutch output member will 
Mewar square wave with hali-cycle duration determined by the pulse width 


modulator input. A Fourier decomposition of this signal yields: 


OO 


jm) t 
Melb) = 3 eee 
Opis = (e0. 
jm, T 
a ) Pe ciggn 4210 
jmo) 
de (Cig). = 
(2 - -1)A ; m=O 
where 82 = clutch output member angular velocity 
@) = switching rate (rad/sec) 
a = imag mitude eiamput imeémber speed 
T = duration of positive half-cycle of PWM output 
ga) SS oyzyrorancigiite baleen 


i —epcmiod oe: — iM output 
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Theodorsen's Function 


Phen ee. 





2 Oe 


imeunemease Of steady state zero positiomerror, T = >: inc 
magnitude of the fundamental component of clutch output velocity will 
4 
be a maximum of Bees : 
al 


Referring again to Fig. 1.2(a) relations between control surface 


output velocity and clutch velocity can be written: 








Zz 
vc. ae (ean is : is 
Q2 n K ik 
f ff 
where Ze = load impedance 
fie 2=) Teauielhon rate 
Se cer pleret requency, 


For switching frequencies above load resonance (25 cps), inertia 


terms dominate and Ze — sJi- 


Mv 
Me COMDINIMg : = = on Il as ais ae 
L f 


J 
Z Ibs 
S (J, + aad 
n 





Pxcmining the relative magnitudes of inertias: 


J 


2? 


3x 10> in lb Bee (Ce Jo Sec tice eoe 


iF 8x 194 in lb eee 


L 


i" 


sance 2a reduction ratio of at least 15:1 will be required, the effects 


of reflected load inertia can be ignored with the result that 


a 


2 J 


njlts K | 





indicating that switching frequency harmonics will be attenuated 


inversely as the square of the order. Coupling this result with the — 





yee 
Chaweebetmer lie ClutCh output spectrum, the foil motion can be 
considered a pure sinusoid at fundamental switching frequency. 


The source SPL can now be calculated 


ie dL, is 
e 7 ( dt 
source 16 ao ‘ 2 
O 
where L= Lo... sin ot = C(k) Ky [o,, sin ot] 
Ky = constant of lift-incidence proportionality 
Ke = 1 pS Oe (MKS units) 


Jl 
C(k) = Theodorsen's function 
Bimee the Seund pressure is a result of vertex shedding, 2 
measure of spanwise ae correlation is necessary. A good measure 
Siepnase correlation is obtained by a strip theory appreach. The foil 
was arbitrarily subdivided into four spanwise sections. The Theodorsen 
function was found to be nearly invariant across the span so that the 
sound field was assumed 100 percent phase correlated. The tip vortex 
baffle increases the reliability of this assumption. 


Comverting (Me mean square pressure to decibel ratio 


SPL ource = «0 logyg (2, ) + P, dB 
ee 
eee 
lees a 
f 
P, contains Ko meCircmeerererence pressure (0,002 ub), (0.707)*, 


16 a@ and oe : 





yy) 


Q2 
If we further assume that — = 100 eM (specified maximum 
CantrolesWmtace specd), then 
SPL = P, - 20 lo es: na 
source. = 2 SOM =) f 


This relationship is plotted in Fig. 3.2 for various ratios of 
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switching frequency to F 
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Fig. 3.2 Source SPL vs Switching Frequency 
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CLUTCH POWER DISSIPATION 


As mentioned in Section 1.4 power is dissipated in the interval 
between the switching instant and the point at which clutch input and 
output members become locked. The instantaneous dissipation is the 
peecduct of pullout torque and slip speed. The time average power 
dissipated will be approximately twice the product of clutch output 
inertia, input member speed squared and switching rate. 

In order to minimize power dissipation, the absolute value of slip 
Melociby must be minimized, Assuming J. ; Jy and switching frequency 


meced (see big, ].2{(s)), the effect om slip velocity of varying K, was 


if 


observed on an analog simulation. Since J, SAUL [sys ene et= inte) 1akive 

binges greater than Jo ; Jy was assumed infinite for purposes of 
Simulation. The system simulated and computing diagram are indicated 
Poeeieure 4, i(a). 

The simulation was carried out on the Philbrick Analog Computer 
pme Wiat, |. Engineering Pmejects Labomatonys Results ameundicated 
fmebieune 4,1(b). The input member velocity is the square wave. 
femeeoweput miemiber Velocity is the initially ramp-hke square wave 
Which tracks the input. The third trace is the amplified difference 
between the two. A small biasing error is evident in scme of the 
Semputations, This is attributed to amplifier drift since if was necessary 
to saturate one amplifier in the loop. 

From the simulation, it may be concluded that the spring Ky 


should be made as stiff as possible without exceeding clutch pullout 


=A - 
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SYSTEM SIMULATED 





COMPUTING DIAGRAM 


Fig. 4.1 (a) Analog Simulation of Clutch Dissipation Phenomena 
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; up to pullout torque 
Fig. 4.1(b) Response of Second Order Clutch Output System 





7 
torque in spring windup. At first glance, no differences may 


be apparent in Fig. 4.1(b). A closer examination reveals that response 
w w 

ime tones = ll cue naboUimnalt chavetor 
sw sw 


= =) 0.478. 








Input member speed and switching rate must be accounted for 


in the final selection of Ke 
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